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A model of an anti-phase modulated d-wave superconductor has been proposed to describe the 
decoupling between Cu-0 planes in 1/8 doped La2-xBaxCu04. Unlike a uniform d-wave supercon- 
ductor, this model exhibits an extended Fermi surface. Within Bogoliubov-de Gennes theory, we 
study the mixed state of this model and compare it to the case of a uniform d-wave superconductor. 
We find a periodic structure of the low-energy density of states, with a period that is proportional 
to B, corresponding to Landau levels that are a coherent mixture of particles and holes. These 
results are also discussed in the context of experiments which observe quantum oscillations in the 
cuprates, and are compared to those for models in which the Fermi surface is reconstructed due to 
translational symmetry breaking in the non-superconducting state and to a model of a Fermi-arc 
metal. 



I. INTRODUCTION 

High temperature superconductivity in the under- 
doped cuprates condenses from a so-called pseudogap 
phase whose properties are distinctly different from those 
of a conventional metal. Below a relatively high tem- 
perature, T*, a gap which may or may not be connected 
to superconductivity starts to develop in the excitation 
spectrum and affects the temperature dependence of all 
transport processes. At lower temperatures, above but 
closer to Tc, superconducting fluctuations in the form of 
a disordered vortex- ant ivortex liquid grow up until long- 
range d-wave superconducting order appears at 

The nature of the pseudogap phase has been the sub- 
ject of much study and debate. One can characterize the 
behavior in this phase as "non- Fermi-liquid-like" , which 
typically means that the sharp fermionic excitations of 
a Fermi liquid are broadened even close to the Fermi 
surface (FS), although the situation in the pseudogap is 
somewhat more complicated. ARPES observes "Fermi 
arcs" - sections of 2D FS which appear to terminate at 
gaps and which become shorter, possibly tending toward 
nodal points, as T is lowered. No sharp quasiparticles 
are observed near the anti-nodal points at (±7r, 0) and 
(0,±^). 

The origin of this non- Fermi- liquid- like behavior is 
hotly disputed. Some studies connect it to resonating va- 
lence bonds or preformed pairs, ^ while others associate 
it with exotic forms of fluctuating or static spatial or- 
der such as charge or spin density wave^^i^ or singlet or 
triplet D-density waves. ^ 

This unsatisfactory state of understanding was com- 
pounded in 2007 by the observation of quantum oscilla- 
tions, first in the Hall resistivity and shortly afterward 
in Shubnikov-de Haas and in de Haas-van Alphen mea- 
surements, below Tc at fields that are comparable to, al- 
though typically smaller than, i^c2-^^^^ Since such fields 
are still weak relative to Fermi energy scales, one might 
expect that, once superconductivity is destroyed by a 
large magnetic field, the underlying resistive state would 
be the same pseudogap state as exists above Tc, that is, 
a non-Fermi liquid. However quantum oscillations are 



normally associated with sharp, closed Fermi surfaces of 
a Fermi liquid. Furthermore, it is not straightforward 
to connect the FS areas determined by quantum oscilla- 
tions with the Fermi arc observed in ARPES. However 
it has been noted that commensurate static translational 
symmetry breaking, due to charge or spin density waves, 
could reorganize the large hole FS of the undistorted lat- 
tice into a number of smaller hole and electron pockets 
and that the small electron pockets could account for the 
quantum oscillations, while sections of the larger hole 
pockets coincide with the Fermi arc— The explana- 
tion for the arcs is then that the spectral weight due 
to the periodic perturbation of the charge-density wave 
(CDW) or spin-density wave (SDW) is large on the arcs 
that are observed by ARPES and small on the remainder 
of the FS hole pockets. Such a result is consistent with 
a simple picture of zone folding due to a weak periodic 
superlattice potential. 

In this paper, we consider a variation of this picture in 
which the periodic superlattice arises from a modulation 
of the d-wave superconducting gap function^^^Ti^i Such 
a modulation has been invoked to explain an interesting 
phenomenon called the 1/8 anomaly which is observed 
in some of the lanthanum cuprates. Most lanthanum 
cuprates exhibit singular behavior in the doping depen- 
dence of various low-temperature properties around 1/8 
doping, which is known to coincide with a charge stripe 
structure with the periodicity of 4 lattice constants. The 
superconducting condensate for the proposed model oc- 
curs at a nonzero wave vector, corresponding to a period 
twice that of a charge stripe structure. ^^^^^ Furthermore, 
the stacking arrangement assumed for this model results 
in a zero Josephson coupling between nearest cuprate lay- 
ers which explains the apparent dynamical decoupling of 
cuprate layers observed in transport measurements of 1 /8 
doped La2-ccBaa;Cu04.^^ In this paper, we study such 
a modulated d-wave gap in the presence of large mag- 
netic fields using lattice Bogoliubov-de Gennes (BdG) 
theory^^ to determine if quantum oscillations, associated 
with Landau level formation, occur. 

The remainder of this paper is organized as follows. In 
Sec. 2, we describe the model for a 7r-striped supercon- 
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ductor. This section shows the density of states (DOS), 
the spectral functions and the FS in zero field for differ- 
ent amplitude gaps for this model. In Sec. 3, we establish 
the generation of Landau levels in the DOS by a magnetic 
field and the effect of doping and of the gap amplitude on 
the Landau level spectra. In Sec. 4, the specific heat is 
calculated to make some connections to the experiments. 
Finally, in Sec. 5, we discuss the results and compare 
them to other models and to relevant experiments. 



II. THE MODEL AND METHOD 

The two-dimensional tight-binding model for a tt- 
striped superconductor is described by the following 
mean-field Hamiltoniar^ 



H = Hq + ^ A{cos(g^x)[4^^^4+i^^^ - 4,^^4+1^^^] 

-cos((7^(x- l/2))[4 .4 



(1) 

] + H.C.} 



where cj.^^ creates an electron with spin cr on site (x, y). 
By setting = 7r/4, the Hamiltonian describes a system 
with a d- wave- type order parameter that has a sinusoidal 
modulation with 8-site periodicity in the x direction, i^o, 
the kinetic part of the Hamiltonian, has the dispersion 
eo = — 2t(cos(A:aj) + cos(/c^)) — At2 (^os{kx) Q^os(ky) — /i in 
k space where t and ^2 are the first and second nearest 
neighbor hopping terms. Due to the periodic modulation 
of the order parameter, the superconducting condensate 
occurs at a non-zero and a particle with wave vector k 
is coupled to ones with wave vectors —k^Qx. We shall 
see that this property of a striped superconductor has 
crucial effects on its low-energy properties. 

There are two possible stable configurations for the 
order parameter of the 7r-striped superconducting model. 
One configuration is the site-centered configuration in 
which the node of the modulation lies on a site. The 
other one is the bond-centered configuration in which the 
node lies on a bond. The calculations in this work are 
done for the latter configuration which is shown in Fig. 
[TJ However, the qualitative behavior of the system in the 
presence of a magnetic field is similar for the site-centered 
case. 

We solve the model by diagonahzing the BdG 
Hamiltonian.^^ The local density of states is defined as 

D{i, E) = 2Y^ [\un(i)\^S{E - En) + \vn(i)\^S{E + En)] 

(2) 

where is the number of sites in one magnetic unit cell, 
N is the number of magnetic unit cells and {un{i),Vn{i)) 
is the position eigenstate of the n-th positive-energy 
state. Due to translational symmetry, the local density 
of the system is the same for all unit cells. The spectral 
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FIG. 1: Position dependence of the pairing gap for the bond- 
centered configuration using color coding on bonds. The cir- 
cles in the middle of the plaquettes specify the positions of 
vortices for a / = 8 magnetic field unit cell whose boundary is 
shown by the dashed line. In the singular gauge, the vortices 
at white (dark) circles are only seen by particles (holes). The 
lower part of the figure shows the varying gap amplitude as a 
function of x. 
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FIG. 2: DOS of a 7r-striped superconductor for various values 
of the pairing gap amplitude. The second nearest neighbor 
hopping for this DOS calculation is set to zero and the chem- 
ical potential /i is adjusted to yield 1/8 doping here and in 
the following figures unless another value is explicitly stated. 
Note the finite DOS at zero energy and the complex structure 
which arises from band folding associated with the strength 
of the periodic interaction as discussed in the text. 



weight function in the extended Brillouin zone (BZ) is 
defined as 



Nl 

A{k, E) = ^J2 {\un{k)f5{E - En) + K{k)\^5{E + i?„) 

n=l 



(3) 
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FIG. 3: The spectral weight (left) and FS (right) for four 
values of the pairing gap A a) 0.05, b) 0.1, c) 0.2 and d) 0.4 
in half of the extended BZ. The colorbar applies only to the 
spectral weight. 



where iXn(^) ('^n(^)) is the Fourier transform of Un{i) 
{vn{i))- In the reduced BZ scheme, one sums A{k^E) 
over the eight coupled k in the extended BZ that can be 
folded back to one point in the reduced BZ. The DOS as 
a function of energy can be obtained from the position 
average of the local density of states or the wave vector 
average of the spectral weight function. 

The DOS of a homogeneous d-wave superconductor 
vanishes linearly at low-energy. In contrast, a tt- striped 
superconductor has a non-zero DOS at zero energy.— The 
low-energy dependence of the DOS for various values of 
the pairing gap amplitude. A, is shown in Fig. [2j For 
simplicity, the second nearest neighbor hopping in the ki- 
netic part of the Hamiltonian has been set to zero and all 
the energy quantities are written in units of the nearest 
neighbor hopping t which is set to 1. The chemical po- 
tential, /i, is adjusted to yield 1/8 doping unless another 
value is explicitly stated. 

It is useful to compare and contrast the A dependence 
of the DOS shown in Fig. [2] to that of the one-electron 
spectral weight. For small A, small gaps open in the 
unperturbed FS segments that can be connected by Qx 
as shown in Fig. [3l^a). Consequently the DOS at zero 
energy does not change significantly with respect to the 
unperturbed case. For intermediate values of A, in the 
approximate range (0.07 < A < 0.13), the pairing gap 
has become strong enough to create two sets of closed 
loops which can be seen in Fig. [3fb) for A = 0.1. The 
Fermi velocity associated with these loops is small and 
consequently they contribute considerably to the DOS 
at zero energy. This is why there is a peak in the DOS 
for A = 0.1. For larger values of A, in the approximate 
range (0.14 < A < 0.25), the loops are gapped out and 
the peak disappears. In this range, the spectral weight 
exhibits Fermi arcs with two small gaps. For even larger 
A, the gaps within the Fermi- arc-shaped spectral weight 
become larger and the shape of the FS in the repeated 
BZ scheme appears as figure-8-shaped loops as shown in 
Fig. ^d) for A = 0.4. 



III. RESULTS IN A MAGNETIC FIELD 

The above results are all for zero magnetic field. A 
magnetic field is incorporated into the model using the 
Franz- Tesanovic singular gauge transformation.^^^^^ In 
this approach, the magnetic unit cell has linear size 
where / is measured in units of the lattice constant. Each 
unit cell has two vortices; one is seen only by particles 
and the other seen only by holes. We position the vor- 
tices at the nodes of the order parameter, as shown in 
Fig. [1] for the case of / = 8. We take / to be an integer 
multiple of 8 which is the period of the order parameter. 
The magnetic field and / are related hy B = (l)o/P where 
(j)o is the flux quantum. For example, taking the lattice 
constant a = 3.8A, / = 32 corresponds to 5 = 28 T. 
Hereafter, we express the magnetic field in terms of 1. 

In this section, we investigate how the DOS structure 
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FIG. 4: Low-energy DOS of a vr-striped superconductor with 
A = O.Olt and /x = —0.226 in the presence of magnetic fields 
of / = 24 (top) and / = 32 (bottom). 



of the model changes as a function of the pairing am- 
plitude A in the presence of a magnetic field. First, we 
consider small values of A where one can expect to un- 
derstand the effect of the interaction based on a simple 
perturbative picture. For A = 0.01, the spectral weight 
in the absence of a magnetic field exhibits only small gaps 
at four points of the unperturbed FS. It is similar to Fig. 
[3fa) except that the gaps are smaller. 

The low-energy DOS structures for A = 0.01 and two 
magnetic fields of / = 32 and / = 24 are shown in Fig. HJ 
The most striking feature of this figure is the appearance 
of Landau levels that are equally spaced in energy with 
the spacing proportional to B. Furthermore, the presence 
of a small perturbative interaction. A, causes the low- 
energy Landau levels to be slightly broadened and also 
partially reflected to the other side of the Fermi energy 
due to particle- hole scattering. In fact, each Landau level 
for A = is split into two peaks with the second peak 
having much smaller weight for small A, as seen in Fig. 
m The sum of the number of states in these two peaks 
equals the degeneracy of a Landau level. 

From a semiclassical point of view, particles can keep 
undergoing Larmor precession by tunnelling through the 
gaps since the gaps are small for A = 0.01. This is the 
so-called magnetic breakdown phenomenon. A particle 
can also be Andreev scattered as a hole into a state of 
—k^q. This process explains the reflected part of each 
Landau level with smaller weight in Fig. HJ This pic- 
ture is motivated by the work of Pippard,^^ who studied 
the cyclotron motion of nearly free electrons in the pres- 
ence of a weak periodic potential that induces gaps in 
the Fermi surface. For this case, when the periodic po- 



tential is weak, electrons can tunnel through the gaps, 
following the unperturbed FS trajectory, or they may be 
Bragg scattered onto a different cyclotron orbit leading 
to broadening. The main difference between Pippard's 
model and the 7r-striped superconducting model is that 
the superconducting periodic potential scatters electrons 
with wave vector k into holes with wave vector —k^q and 
vice versa. Thus electrons either tunnel through gaps in- 
duced by the periodic potential or scatter into hole states. 
Note that magnetic breakdown occurs even if the magni- 
tude of the gap in the FS is larger than hujc^ 

For intermediate values of A (0.07 < A < 0.13), for 
which the FS has well-separated segments (see Fig. [3l^b)), 
we do not observe clearly defined Landau levels. This 
may be the result of broadening and the close spacing of 
Landau levels due to the large density of states. Further- 
more, multiple Fermi surfaces may each give rise to their 
own sets of Landau levels which are unresolved. 




FIG. 5: DOS for A = 0.25 and magnetic field of / = 32 shown 
as a function of positive and negative energies separately. The 
band structure spans energies from — 4 — /x to 4 — /i. However, 



the DOS is only shown in the 


-1 < E < 1 range. 
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FIG. 6: Low-energy DOS for A = 0.25 and magnetic fields of 
/ = 40 (top) and / = 32 (bottom). 

For larger values of A (0.14 < A < 0.3), the shape of 
the FS is simpler. In this range, the spectral function has 
significant weight on the parts of the FS that resemble 
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the Fermi arcs observed in ARPES experiments. The 
DOS for A = 0.25 and / = 32 is shown in Fig. 6 for 
positive and negative energies up to E = 1 separately. 
Remarkably, we again observe periodic behavior of the 
low-energy DOS as a function of E with a spacing that 
varies linearly with B. This is illustrated in Fig. [6] for 
A = 0.25 and two values of the magnetic field, / = 40 
and / = 32. Note the splitting of each Landau level into 
a strong and weak peak seen for small A, Fig. IH does 
not occur in this larger A range, where the original large 
FS is not accessible to the quasi-particles. 

In Fig. [6l the DOS has a minimum, or possibly a very 
small gap, at £^ = for / = 40. However, for / = 32 
it appears that the two Landau levels closest to £^ = 
are joined together and the DOS at £^ = has a nonzero 
value. In general, we find that, for / = 8m where m is 
an integer, if m is even, the DOS at £^ = is nonzero 
and if m is odd, the DOS is zero at £^ = 0. This is a 
commensurability effect that results in oscillation of the 
DOS at £^ = as a function of / or 1/ ^/B and is discussed 
further in appendix [Al 
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FIG. 7: Low energy Landau level spacing as a function of 
l/Z^ for A = 0.25 and /i = -0.3. The spacing is defined as 
E(N)/N where E{N) is the minimum in the DOS between 
the N-th and N + 1-th Landau levels and is shown for N — 2 
(triangle) and N — 1^ (circle). The line is a linear fit to the 
data that goes through the origin. 

We have calculated the spacing of the low energy Lan- 
dau levels for a wide range of fields for A = 0.25 and 
/i = —0.3 as shown in Fig. [71 The Landau level spac- 
ing can be defined as E{N)/N where E{N) is the mini- 
mum in the DOS between the A^-th and N + 1-th Lan- 
dau levels and is essentially independent of N provided 
E{N) < 0.5A. Fig. [71 demonstrates the spacing as a 
function of 1/P oc B for N = 2 and N = 10. The slope 
of the Landau level spacing versus B is inversely propor- 
tional to the DOS dit E = 0. By comparison, we find that 
the slope is about half as large and the DOS at £^ = 
about twice as large for A = 0. 

Furthermore, the number of states in each peak is 
nearly the same as that of a Landau level. In general, 
in the presence of a magnetic field, the n-th peak on the 
left of £^ = can have a degeneracy slightly different 



from a Landau level degeneracy. However, the n-th peak 
on the right compensates so that the number of states of 
the two peaks together is always twice that of a Landau 
level. This shows that the Landau levels are a coherent 
mixture of particles and holes together and the particle- 
hole scattering is playing a role in the formation of the 
Landau levels. The reason for the small difference of the 
number of states in each peak from the exact degener- 
acy of a Landau level is that the low-energy DOS in the 
absence of a magnetic field is asymmetric around E = 
except at half filling, as shown in Fig. [51 




FIG. 8: The DOS structure in the absence of a magnetic field 
for A = 0.25 and two dopings. Note the asymmetry at low E 
for 1/8 doping. 

It is worth contrasting the behavior of the 7r-striped 
superconductor. Fig. [6l to the DOS structure of a homo- 
geneous d-wave superconductor. For the latter at half- 
filling. Landau levels are formed in the low-energy DOS, 
but the Landau levels are not equally spaced and the 
spacing scales as ^/B around E = 0.^^ This is a conse- 
quence of the nodal behavior at the Fermi energy. There- 
fore, quantum oscillations periodic in 1/B are not ex- 
pected for a d-wave superconductor. For the remainder 
of this paper we refer to each peak of the type shown 
in Fig. [6l (that is, equally spaced with a spacing propor- 
tional to B) as a Landau level. The fact that there is only 
one set of Landau levels and the number of states in each 
peak is equal to that of a Landau level indicates that all 
parts of the FS participate in the formation of low-energy 
Landau levels. As noted above, the case of large values 
of A, Fig. [6l is different from the case of very small val- 
ues of A, Fig. [H in which the sum of the degeneracy of 
the two peaks equals the degeneracy of only one Landau 
level. 

Although Landau levels are observed in the low-energy 
DOS in the large A regime (0.14 < A < 0.3), for very 
large A (A > 0.35), the Landau levels disappear. Fig. [9l 
shows a comparison of the low-energy DOS for A = 0.2 
and A = 0.4 in the presence of a magnetic field / = 48. 
The DOS for A = 0.4 shows no Landau levels and a sharp 
peak at ^ = which can be attributed to the commen- 
surability effect. This is further discussed in appendix 
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FIG. 9: Low-energy DOS for / = 48 and A = 0.2 and A = 0.4. 
Landau levels are suppressed for A = 0.4 but a sharp peak 
around ^ = appears. 

El 

As A is increased, large gaps appear within the Fermi 
arcs. This is seen in the spectral function, shown in Fig. 
[31(d) for A = 0.4 and zero magnetic field. Assuming that 
the semiclassical process of the formation of the Landau 
levels^^ involves tunnelling of particles (holes) across the 
gaps within the Fermi arcs, one expects the magnetic 
breakdown phenomenon not to occur if the gaps are too 
large. This may explain why the Landau levels are sup- 
pressed for very large A. 




FIG. 10: Density of electrons versus — /i for the magnetic field 
of / = 16 and two cases of A = 0.25 and A = 0. Unlike A = 0, 
the density does not exhibit a stepped behavior for A = 0.25. 

Landau-type quantum oscillatory behavior has previ- 
ously been discussed in the context of a particular model 
of a 'Fermi- arc metal'. In that model, parts of the FS of 
a metal are artificially gapped out by restricting super- 
conducting pairing to the antinodal regions of momen- 
tum space in order to get a FS that consists of Fermi 
arcs. The 7r-striped model, which is based on a specific 
microscopic mechanism and has no such restriction, dif- 
fers from the Fermi-arc metal of Ref. 33 in that the 
one-electron spectral function has non-zero (but possibly 
very small) weight along continuous lines in /c-space. 

At a more basic level, the behavior of a 7r-striped su- 
perconductor is strikingly different from that of a metal, 
in spite of the fact that both exhibit a Fermi surface. In a 
metal, the particle density n versus fi exhibits a stepped 



behavior in the presence of a constant magnetic field. 
In contrast, the particle density in a 7r-striped supercon- 
ductor changes smoothly as a function of /i as shown 
in Fig. [To] for a large magnetic field of / = 16. Fur- 
thermore, we find that the low-energy DOS behavior of 
the TT-striped superconductor is rather insensitive to the 
change in ja. In other words, no oscillatory behavior of 
the DOS at £^ = is observed as fi is varied except for 
finite size effects. This is in contrast to the result for 
the simple Fermi-arc metal model. However, quantum 
oscillations are induced by changing the magnetic field, 
not the chemical potential, and consequently could still 
be observed for a 7r-striped superconductor. 

IV. SPECIFIC HEAT 

In this section, we provide specific heat calculations 
in the absence and presence of a magnetic field to make 
connections to experiments on the cuprates. In advance, 
we note that the dependence of the Sommerfeld co- 
efficient, 7, in the cuprates is not present in the 7r-striped 
superconducting model as there is a finite DOS at £^ = 0. 
However, it will be shown that some features of the spe- 
cific heat in the cuprates are consistent with this model. 
One can calculate the specific heat by using the relation- 
ship c = T|^. For a system of quasiparticles, the entropy 
is given by^^ 

S = -kB Y.[U In /p + (1 - fp) ln(l - U)] (4) 

pa 

where a is the spin state and fp is the Fermi-Dirac dis- 
tribution function given by 



and ep is the energy of the quasiparticle associated with 
the state p. Converting the sum over states in Eq. [J] 
to an integral over energy brings in the DOS. It should 
be noted that, in this study, the DOS is not calculated 
self-consistently as there is no microscopic Hamiltonian 
defined. Furthermore, it is assumed that the magnitude 
of the pairing interaction is fairly constant at low temper- 
atures so that the quasiparticle spectrum is unchanged as 
temperature increases. 

The specific heat at zero field for a 7r-striped supercon- 
ductor as a function of temperature is shown in Fig. [11] 
for A = 0.25 at 1/8 doping. The slope associated with 
the linear behavior is 7 = 0.70k%t~^ per site and is di- 
rectly proportional to the DOS at £^ = 0, which is 0.21t~^ 
per site. The slope is about half of that of A = at 1/8 
doping. The low-energy specific heat of various fields as 
a function of temperature is also shown in Fig. [11] for 
A = 0.25 and /i = —0.3. As expected, all the curves con- 
verge to that of zero field as the temperature increases. 
However, at very low temperatures, the specific heat be- 
havior for different fields is significantly affected by the 
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commensurability effect. This is seen in the nearly zero 
slope of the curves for odd m (recall / = 8m) as T ^ 0. 




T(t/kJ X10-' 



FIG. 11: Specific heat in the presence of various fields as 
a function of temperature for A = 0.25 and fi — —0.3. The 
behavior of the curves at very low temperatures is significantly 
affected by the commensurability effect. The heavy line shows 
the specific heat in zero field for A = and fi — —0.225 
corresponding to 1/8 doping. The slope associated with the 
linear behavior is about two times that of the slope for A = 
0.25 in zero field (noted by / = oo). 

Low temperature electronic specific heat measure- 
ments of a cuprate^^ indicate a relatively large DOS at 
E = ^ which can not be explained by the presence of 
disorder in a d-wave superconductor. Taking the lattice 
constant of a typical cuprate to be a = 3.85A, the specific 
heat effective mass becomes m*/m = 0.34eV/t. A rather 
wide range of values has been used for t?^'^^ Within a 
simplified nearest-neighbor hopping only model as used 
here, we obtain m*/m = 1.36 for t ^ 0.25eV.^^ This 
value corresponds to 7 ~ 1.98 mJ- K~^-mol~^ which 
is consistent with the specific heat measurements for 
cuprates in the absence of a magnetic field, 7 ~ 1.85 
mJ- K~^-mol~^J^ 

Riggs et al. have studied the low temperature spe- 
cific heat as a function of magnetic field up to very high 
fields (50 T) and observed quantum oscillations.^^ This 
allowed them both to measure the magnitude of the spe- 
cific heat in what is presumably the normal state, and 
also to determine the cyclotron effective mass associated 
with the quantum oscillations. Then if one assumes the 
normal state has broken translation symmetry, model- 
ing the arrangement of electron and hole pockets in the 
Brillouin zone and using the measured cyclotron effec- 
tive mass, one can estimate what the specific heat should 
be. The result is much larger than the specific heat that 
they observed. ^^^^^ The same problem was also noted in 
a theoretical study^^ based on FS reconstruction where 
the calculated specific heat was larger than the measured 
value. 

In our calculations for a 7r-striped superconductor, it 
was found that, even though there exist several FS pock- 
ets in the BZ, a single set of Landau levels is observed 
above and below the Fermi energy. The relation between 



the slope of Fig. [71 defining the cyclotron effective mass 
obtained from the spacing of Landau levels, and the DOS 
at £^ = is the same as for the A = case. This implies 
that, as for the A = case, the cyclotron effective mass, 
mc, is equal to the specific heat effective mass, m*, for 
large values of A. Consequently, the quantum oscillations 
in the specific heat and the magnitude of the specific heat 
which is observed in Ref. 12 could be consistent with the 
behaviour of a 7r-striped superconductor state induced 
by large magnetic fields, rather than a striped metallic 
state with no pairing gap as is often assumed. However, 
as noted earlier, the ideal 7r-striped model (with no uni- 
form (i-wave component) is not expected to give a 
background, which also appears to be a feature of the 
experiments.^^ In addition, since we cannot study small 
changes in B we make no prediction about the spacing 
of the observed quantum oscillations. 



V. DISCUSSION AND CONCLUSION 

In this work we have studied the requirements for hav- 
ing quantum oscillations in a model of a 7r-striped su- 
perconductor. For a large range of values of the pairing 
interaction, the FS corresponds to closed loops while the 
one-particle spectral function exhibits Fermi arcs in k- 
space. Our main finding is that Landau levels are seen in 
the low-energy DOS of the 7r-striped superconductor in a 
large range of the magnetic field, which indicates the pos- 
sibility of quantum oscillations. We find that low-energy 
Landau level formation persists even though particle and 
hole levels are mixed by the pairing interaction. Other 
theoretical studies of quantum oscillations in the cuprates 
are typically based on FS reconstruction of a metallic 
state and involve multiple pockets and frequencies. ^-1^ 
Furthermore, the pockets associated with those studies 
are located where the ARPES experiment shows a large 
pseudogap. By contrast, the 7r-striped superconductor 
exhibits a unique low-energy Landau level set that is only 
due to the Fermi arc part of the spectral weight function. 

Since our numerical studies are restricted to satisfy- 
ing / = 8m, we cannot change the magnitude of the 
magnetic field continuously or in small steps. In addi- 
tion, the Landau levels are located symmetrically around 
£^ = for the discrete values of the magnetic field that 
we can study. As a result it is not possible, from these 
calculations, to find the FS area associated with quan- 
tum oscillations that would be observed by conventional 
experimental methods. However, we can make conjec- 
tures about FS areas that might be observed, based on 
our analysis. We expect that any semiclassical trajectory 
describing the formation of Landau Levels should have 
the following characteristics: 1) The trajectory should 
use all parts of the FS. 2) Andreev scattering needs to 
occur at least at one point during the Larmor precession 
because Landau levels are a coherent mixture of parti- 
cles and holes. 3) Magnetic breakdown is likely involved 
in Landau level formation because once the gaps within 
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the Fermi arcs become too large, the Landau levels dis- 
appear. 

It is also useful to compare the behaviour of the tt- 
striped superconductor to the Fermi-arc metal, in which 
a new mechanism for quantum oscillations is proposed 
that is not based on FS reconstruction.^^ For that model, 
it was shown, based on a semiclassical approach, that the 
frequency of quantum oscillations is proportional to the 
Fermi arc length. In the 7r-striped superconductor, the 
fact that a quasiparticle with wave vector k is coupled 
to ones with wave vectors —k — Qx and —k-\-qx pro- 
vides a different scattering mechanism which changes the 
semiclassical motion of quasiparticles. Consequently, the 
semiclassical trajectories of the two studies are expected 
to be different. Although Landau level formation in the 
Fermi-arc metal resembles what we have seen in the tt- 
striped model, the DOS at £^ = shows an oscillatory 
behavior as a function of /i for the Fermi-arc metal which 
is inconsistent with our study. 

We note that the observation of quantum oscillations 
corresponding to small Fermi surface pockets supports 
the scenario of translational symmetry breaking and 
Fermi surface reconstruction, whether due to charge or 
spin density waves or to modulation of the d-wave gap, 
as discussed in this paper. Indeed recent experimental 
results support the connection between stripe formation 
and quantum oscillations.^^ Our calculations show that 
a modulated d-wave superconductor can support Lan- 
dau levels and quantum oscillations but we are unable to 
make detailed comparisons to quantum oscillation exper- 
iments because of the restriction to commensurate vor- 
tex lattices. One might expect there to be observable 
differences between quantum oscillations in the presence 
of charge or spin stripes and superconducting stripes, 
due to the Andreev reflection and particle-hole mixing 
involved in the formation of Landau levels in the latter 
case. Therefore, it would be of interest to study modu- 
lated superconductivity within a framework which allows 
the magnetic field to be varied continuously to more di- 
rectly connect to the quantum oscillation experiments on 
the cuprates. Possible approaches would be to use ran- 
dom vortex lattices, as was done by Chen and Lee,^^ 
or to develop a semiclassical approximation that allows 
magnetic unit cells of arbitrary aspect ratios. 



Appendix A: Commensurability effects 

In this appendix, we discuss the commensurability ef- 
fect mentioned in Sec. IIIII We found that the DOS at 
E = exhibits a periodic behavior as a function of / or 
\l\fB for large values of A. If m is even, the two vortices 
in a magnetic unit cell are in perfectly equivalent posi- 
tions with respect to the spatially modulated gap and the 
DOS at £^ = is nonzero. In contrast, for odd m, the 
gap on the right of one vortex is positive but on the right 
of the other vortex is negative as can be seen in Fig. [TJ 
The DOS at £: = is zero for odd m. 

The same kind of commensurability effect is also seen 
in a uniform d-wave superconductor.^^ There, due to 
strong internodal scattering, the DOS at zero energy ex- 
hibits a periodic behavior as a function of k^l where k^ 
is the /c-space half distance between the nearest nodes of 
the d-wave superconductor. Specifically, depending on 
whether n is odd or even in k^l = 7rn, the DOS around 
£^ = shows a linear or gapped behavior. For a 7r-striped 
superconductor, the relevant /c-space half distance is tt/S 
which leads to a periodicity of = 16 for the DOS at zero 
energy as a function of L So the commensurability effect 
seen in a 7r-striped superconductor is most likely due to 
interference effects. This suggests that the nonzero DOS 
at £^ = for even m is due to constructive interference 
of particle and hole waves, while the gapped behavior for 
odd m is due to destructive interference. 




5 10 15 20 25 30 35 40 45 



FIG. 12: Local density of electrons due to the low-energy 
states within O.OOlt of ^ = for / = 48 and A = 0.4. 
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For very large A, a sharp peak develops near £^ = 
for even m only, as shown in Fig. [9l It appears that the 
origin of the peak can be traced back to the non-zero DOS 
at ^ = for smaller A, and consequently is related to the 
commensurability effect. The fact that the low-energy 
Landau levels disappear when the peak at £^ = is sharp 
suggests that the commensurability effect is competing 
with the Landau level formation. Fig. [T2l shows the real 
space representation of the the states under the sharp 
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peak at £^ = where a pattern of stripes of low and 
high particle density is clearly visible. On average, the 
ratio of the density of the higher density stripes to the 
lower density is 2.35 for A = 0.4, which is larger than the 
value 0.2 in the absence of a magnetic field (see appendix 
[B]) . It is not clear whether these effects are due to the 
existence of a strong commensurability effect or due to 
the presence of large gaps within the Fermi arcs. 

Appendix B: Periodic Andreev state 




FIG. 13: High (red stripe) and low (blue stripe) density struc- 
ture of the low-energy particles relative to the modulated d- 
wave gap. High (low) density is indicated by a red (blue) 
stripe. 

In this appendix, a type of Andreev state that is seen 
for the low-energy particles (holes) in the absence of a 
magnetic field and persists in the presence of a magnetic 
field is discussed. We have already examined how the 
spectral weight of the low-energy states of a 7r-striped 
superconductor changes as the pairing amplitude A is 
varied. We can also look at the real space representation 
of these states. Our main finding is that, for large enough 



values of A where the shape of Fermi arc is assumed in 
the spectral function (see Fig. [3]), the real space rep- 
resentation of the low-energy states exhibits a periodic 
stripe structure with the periodicity of four lattice sites. 
The stripe structure corresponds to higher and lower den- 
sity of low-energy electrons and holes. Each stripe has 
a width of two lattice constants and the higher density 
stripes are located exactly where the order parameter 
is minimum as shown in Fig. [131 The density ratio of 
the higher density stripes to the lower density ones in- 
creases as A increases. The ratio is approximately 1.5 
for A = 0.2 and 2 for A = 0.4. 

This stripe structure occurs not only at the Fermi en- 
ergy, but also at energies near Fermi energy. The origin of 
these stripes is simple. They are formed due to the con- 
structive interference of the low-energy electron (hole) 
waves with their doubly Andreev scattered counterparts. 
Since the probability of being Andreev scattered twice 
increases with A, the density difference between higher 
and lower density stripes increases accordingly. One can 
think of the stripe structure as a periodic Andreev state. 
For a wave vector at the end-points of the Fermi arcs, 
the doubly scattered wave vector is also located at the 
end-point of another Fermi arc. This means that the 
weights of both interfering waves are large and as a re- 
sult they contribute considerably to the formation of the 
stripe structure. Moving toward the center of the Fermi 
arc, the coupling interaction decreases and the periodic 
Andreev structure is less likely to be formed. This is why, 
as mentioned in Ref. 24, the stripe structure is mainly 
due to states near the end-points of the Fermi arcs. 
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